in 



X 



Thermodynamic properties of spherical nuclei * 

A.M. Kamchatnov 

Institute of Spectroscopy, Russian Academy of Sciences, Troitsk, Moscow Region, 14-2190 Russia 

February 8, 2008 



m 
o 
o 

(N 
o . 

D . Abstract 



The effect of the residual interaction between the nucleons (quasiparticles) on the 
thermodynamic behavior of spherical nuclei is considered. Thermodynamic relations 
are derived for the simplest model of the residual interaction that has the necessary 
macroscopic structure previously established. The range of applicability of the theory 
is discussed. 
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1 Introduction 

m 

The nuclear physics data accumulated up to now indicate that there is an essential difference 
between spherical and nonspherical nuclei. This difference does not reduce merely to "nuclear 
shape effects" as such, but intimately involves the internal structure of the nucleus. The 
shell oscillations of nuclear masses [1] can serve as an example of the phenomena associated 
with the difference between the structures of spherical and nonspherical nuclei. The shell 
oscillations are due to the ordered grouping of one-quasiparticle levels in spherical nuclei that 
arises because of the existence of the orbital angular momentum quantum number I of the 
quasiparticle [1,2]. There is no such ordered grouping in nonspherical nuclei, and the trend 
of the mass of such nuclei as a function of the number of nucleons is entirely different. This 
difference is manifested on a plot of nuclear masses versus number of nucleons as a break 
in the curve corresponding to a phase transformation from the normal phase (nonspherical 
nuclei) to the "magic phase" (spherical nuclei) [3]. 

The difference between the internal structures of spherical and nonspherical nuclei does 
not affect only the ground states: the thermodynamic properties of the two phases are also 
different. This difference is manifested, for example, in the dependence of the form of the 
spectra of the 7-ray cascades emitted by heated nuclei on phase of the compound nucleus 
[4,5]. The energy spectrum of a nonspherical nucleus is not distorted by the residual in- 
teraction and is essentially of the Fermi-liquid type. Spherical nuclei, on the other hand, 
apparently cannot exist without the residual interaction (this was proved in [6] for a very 
simple nuclear model), and this interaction has an important effect on their thermodynamic 
behavior. As was explained in [2], the residual interaction must have a certain macroscopic 
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structure, the only interactions between the quasiparticles that are compatible with the ob- 
served pattern of the shell oscillations of nuclear masses being those that lead to a quadratic 
dependence of the width of the diffuse region of the Fermi distribution on the orbital angular 
momentum I of the quasiparticle: 

5e oc l 2 (Z + 1/2) 2 . (1) 

In this paper we consider the thermodynamics of spherical nuclei on the basis of the 
simplest model of the residual interaction that has the required macroscopic structure. 



2 The Racah-Mottelson model 

The Racah-Mottelson model was used in [2] to describe the residual interaction: the Hamil- 
tonian for the interaction between quasiparticles on the same j level is given by 

H L = ~ G 3 Yl a ln' a - m ' a - m a m , (2) 

m,m'>0 

where a) m and CL'fji £1X6 the creation and annihilation operators for a quasiparticle whose 
angular momentum has the z projection m. The eigenvalues of this Hamiltonian (see, e.g., 
[7]) are 

EL = -GM^-b,-s 3 + l), (3) 

where 2Qj — 2j + 1, bj is the number of interacting pairs on level j, and Sj is the seniority 
quantum number. 

To construct the thermodynamics of spherical nuclei, we need a formula for the degree 
of degeneracy T(s, Q) (we shall drop the subscript j in the rest of this section) of the Racah- 
Mottelson energy levels (3). The symmetry of the zero-seniority ground-state wave function 
(see, e.g., [8]) is such that the wave function can be constructed from single-particle wave 
functions in only one way. Hence the s = ground state is not degenerate: r(0, fl) = 1. One 
more quasiparticle can be added to the s = state in ways (Cf n = (2f2)!/s!(2fi — s)\ 
is the number of combinations of 2Q things taken s at a time), and we thereby obtain an 
s = 1 state, so that T(l, Q) = C\ n . Two quasiparticles can be added to the s = state with 
2b quasiparticles in C| n ways, but among them there is one that describes the ground state 
of 2(b + 1) quasiparticles, so that T(2, Q) = C| n — C® n . Arguing further in the same way we 
can prove the following general formula for the degree of degeneracy of the Racah-Mottelson 
energy levels by mathematical induction: 

T(s,n) = c s 2n -c s ^ 2 , o< s <n. (4) 

For the case O > 1, s > 1, we can use Stirling's formula n\ = ^/2^^nexp(n\nn — n) to 
rewrite Eq. (4) in the form 



r(*, n) = 4 ( ^_ ^ J irs{ m_ s) e Mm in(2n) - 8 in 8 - (2Q - 8 ) in(2n - *)] . (5) 
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Now let us find the number T(M, s, ft) of states with seniority s and angular momentum 
projection M. It is clear that the sum of T(M, s, ft) over all possible values of M is T(s, ft). 
In order to find a way of evaluating r(M, s, ft) let us consider the matrix of Hamiltonian 
(2) in more detail. The dimension of this matrix is the number of ways in which TV 
fermions can be distributed among 2ft cells. (It is easily seen that the sum of T(s, ft) over 
all seniorities s possible for the given N is C^.) The Hamiltonian matrix for the C^q basis 
functions breaks up into blocks corresponding to different values of the quantum number 
s. However, the block for the states with a given seniority breaks up in turn into still 
smaller blocks, since the Hamiltonian conserves the angular momentum projection M, so 
that different values correspond to different blocks. Finally, the block for states with given 
fixed values of s and M also consists of smaller blocks, since different sets of unpaired particles 
can give rise to the same angular momentum projection M. (Here the term "pairing" is used 
in a purely kinematic sense, namely, as an expression of the fact that two quasiparticles lie 
in cells corresponding to the same absolute value but opposite signs of the z projection of 
the angular momentum.) For definiteness let is consider the case of even s (and hence even 
N). Then the "smallest" blocks in the matrix for Hamiltonian (2) describe states in which 
there occur some number 2k of unpaired quasiparticles giving a total angular momentum 
projection M, and N — 2k paired quasiparticles, which have M = and seniority s — 2k. 
The number of such blocks is equal to the number C(M, 2k, ft) of ways in which the angular 
momentum projection M can be expressed as the sum of 2k single-particle projections taken 
from the set { — j, + — 1, j} (no two single-particle projections equal in modulus 

but opposite in sign occurring in the sum), while the dimension of each block is equal to the 
degree of degeneracy T (s — 2k, ft — 2/c) of a state of paired quasiparticles that are distributed 
among ft — 2/c free cells and have seniority s — 2k. Thus, the sum of the dimensions of all 
the blocks for states with angular momentum projection M and seniority s is 

s/2 

r(M, s, ft) = °( M > 2k > fi ) r o(s -2k,n- 2k). (6) 

fc=0 

A similar expression can also be written for odd s. 

Direct diagonalization of the Hamiltonian matrix for small values of ft and iV leads to 
values of T (s, ft) that can be expressed by the formula 

r ( s ,n) = cf-cf-\ (7) 

By definition, To(s,ft) is meaningful only for even s. If we assume that formula (7) is valid 
for all even seniorities < s, we can prove that it is valid for seniority s + 2 in much the same 
way as Eq. (4) was proved: hence it follows by mathematical induction that (7) is valid for 
all even seniorities < s < ft. 

Now let us seek the number C(M, s,n) of ways in which s unpaired particles can give a 
total angular momentum projection of M. Since the formula for T(M, s, ft) will be used later 
to describe macroscopic phenomena for which it is essential that j ^> 1, we can set j = ft 
(also see [2]), so that the single-particle angular momentum projections are the elements of 
the set {—ft, —ft + 1, . . . , — 1, 1, . . . , ft — 1, ft}. Now we shall construct a generating function 
for the unknown numbers C(M.s, ft). To do this we consider the sum 

E(0 £i (c 2 r...(O en (8) 
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taken over the values 



• •; e n = -1,0,1 (9) 

N + ... + M = s. (io) 

The coefficient of ( M in the sum (8) is C(M, s, Q). To take condition (10) into account, we 
make use of the identity 



ei = -1,0,1; e 2 = -1,0,1; 
under the condition that 



l 

2tH 



dz 



|ei| + |£2|+...+|en|-s+l 



Z 
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l £ n| 



s, 



:n) 



in which the integration contour goes around the origin. If we multiply (11) by (8), we can 
sum over all the values (9) independently, without regard for condition (10). As a result we 
obtain the following expression for the generating function: 



M 



1 

2ni 



dz 



n=l 



n i+ 



;i2) 



If we make the substitution ( = e lt in (12), the sum over M becomes a finite Fourier series. 
Using the orthogonality relation 



(13) 



cos Mt cos Lt dt = 7iSml, 
we obtain the following expression for the coefficients: 



1 f dz f n r-l ( 2 

C(M,s,tt) = <b —z s / cosMt Yl ( 1 + -cosnt 

(2ti) i J z J_ w n=i \ z 



(14) 



Now we expand the product of the expressions in parentheses in powers of t and extend the 
integration over the entire real axis, obtaining: 



c(M,s,n) 



(2vr) 2 V fi 3 



37r / dz i 

— (p —Vz + 2 



z s exp 









exp 







3(^ + 2) m2 



4fi 3 



(15) 



This operation is legitimate since the integral over t converges for t ~ ^ (z + 2)/fi 3 and the 
characteristic values of z (see below) are such that z + 2 = 2Q/s, so that t ~ 1/fly/s. We 
evaluate the remaining integral over z by the saddle-point method, choosing the integration 
contour so that the convergence of the integral will be due mainly to the function 

exp[/(t)] = z s exp[ttln(l + 2/z)] = exp[(s - Q) kiz + Qln(z + 2)]. (16) 

At the saddle-point we have 

2(Q - s) d 2 f 



zo = 



dz 2 



ZD 



4Q(Q-s)' 



(17) 
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and we must therefore integrate over the contour z = z + iy. As a result we obtain 



3 1 / O 

c(M, s ,n)= ' ■ 



27rfi^s V 27rs(fi - s) (18) 
x exp[s In 2 + Q In ft - s In s - (ft - s) ln(ft - s)] exp(-3M 2 /2sft 2 ). 

It is easy to show that the condition for the validity of Eq. (18) is that s ^> 1. 

To calculate T(M, s, Q) it remains to substitute expressions (7) and (18) into (6) and 
perform the summation. Approximating the right-hand side of (7) with the aid of Stirling's 
formula and replacing the summation by integration, we obtain 



= /A 3 ^ [ a/2 exp[<K*Q]exp(-3MV4fcn 2 )dfc 
' ' V27T 2VU Jo k{Q - s/2 - k)y/(s/2 - k){Q - s/2 - k) ' 

where 

(j>(k) = tt\nQ-2k\nk- (s/2 - k) ln(s/2 — k) — (Q, — s/2 - k) ln(fi - s/2 - k). (20) 
The function <p{k) has a maximum at the point 

k = s/2 - s 2 /4Q, (21) 

and at that point 

(d 2 <j)/dk 2 )\ ko = -I6{l 3 /s 2 (2tt - s) 2 . (22) 

Thus, if 

s 2 (2n - sf/n 3 < s 4 /n 2 , (23) 

then cf)(k) will have a sharp maximum within the integration range, so we can replace the 
other slowly varying functions by their values at k = k and integrate over k from — oo to 
00. As a result we obtain the following final expression: 

r(M, s, Q) = A ^ S \ exp[(2Q) ln(2Q) - sins - (2fl - s) ln(2Q - s)} . . 

7r s(2\l — s) (24) 

x exp[-3M 2 /^s(2fi- s)]. 

Condition (23) for the validity of Eq. (24) is satisfied when s ^> y/Ti. It should be noted, 
however, that the exponential factors in (24), which are essential for the thermodynamic 
applications, are determined only by the position (21) of the maximum, so that the condition 
s ^> is essentially just a condition for the validity of the preexponential factor in (24); 
the exponential factor, on the other hand, is valid for s ^> 1. On integrating (24) over M 
from — oo to oo we again obtain Eq. (5). 

With that we conclude our study of an individual j level of the Racah-Mottelson model 
and turn to the application of the model to the nucleus as a whole. 
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3 The thermodynamics of spherical nuclei 



Let us calculate the function describing the distribution of quasiparticles in a nucleus heated 
to the temperature T. The preexponential factor in (5) can be neglected in calculating the 
entropy of a single j level with the formula Sj = lnT(sj, Qj). Then we obtain the following 
expression for the entropy of the nucleus as a whole: 



5> = -^2^ 



2fL 2a 



+ 1 



2fL 



In 1 



2fL 



(25) 



To find the function describing the distribution of the quasiparticles among the levels, we 
must minimize the free energy 



F = E - TS, E = J2 {(% + s j )e J - GM^j ~ bj ~ 



(26) 



(here Sj is the energy of a quasiparticle on level j without the residual interaction) under 
the condition that the number of quasiparticles, 



(27) 



remains constant. In view of the macroscopic nature of the problem under study we have 
neglected the "1" in Eq. (3) for the interaction energy (see [2]). 

Condition (27) can be easily taken into account by the Lagrange multiplier method. By 
varying the function F' = F — SfN (e/ is the chemical potential) with respect to bj and Sj 
we obtain the set of equations 



w ■ 



1 

2 



£ f 



GjVtj 



w 



<^'P I j- I +1 



(28) 



for distribution functions Wj = bj/Qj and = Sj/2Qj in the region bj + Sj < Qj, bj ^ 
(region II on Fig. 1). On substituting the expression for from the first of Eqs. (28) into 
the second we obtain the following equation for the distribution function u^: 



cxp 



2T y 



(region II). 



(29) 



On increasing ej, we pass from region II to region III, where there are no paired quasiparticles, 
i.e., bj = 0. On taking this additional condition into account in the Sj variation we obtain 
the following expression for the distribution function w*: 



exp 



T 



+ 1 



-l 



(region III), 



(30) 
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i.e., an ordinary Fermi distribution. On decreasing Ej we pass from region II to region I, 
where all the vacancies permitted by the Pauli principle are filled, i.e., bj + Sj = Qj. This 
auxiliary condition gives a Fermi distribution for the holes: 



to, 



exp 



T 



+ 1 



(region I). 



(31) 



The equations for the boundaries between regions I, II, and III are obtained by equating the 
expressions for the distribution functions on the two sides of the boundary. This gives the 
equation 



T IT IT 

for the boundary between regions II and III. We denote the solution of this equation by 



(32) 



T 



7 



2T 



For the boundary between regions I and II we obtain 

£ j ~ £ .f 



T 



= -7 



2T 



(33) 



(34) 



Since the solution of Eq. (29) is independent of Ej, its value within region II is equal to its 
value on the boundary (for constant Qj) or to the values of functions (30) and (31) on the 
corresponding boundaries (33) and (34), so that 



cxp 



7 



2T 



+ 1 



(region II). 



A graph of the function 7(2;), defined implicitly by the equation 

7 = xtanh(7/2), 



(35) 



(36) 



is shown in Fig. 2. We have 7(2;) = ^/6(rr — 2) as x — > 2 + 0, and 7(2;) = x as x — > 00. 

Equations (30)- (36) completely determine the distribution function Wj and make it pos- 
sible to calculate the entropy of a spherical nucleus from Eq. (25). 

It is convenient to do the calculation in the variables p and (3 introduced in [6,2] (p = kR, 
k is the wave number of the quasiparticle, and R is the radius of the nucleus): 



p[co S/ 3 - (tt/2 - 0) sin (5\ = n(n + 3/4), 

P — arcsin(Z/p), 1 = 1 + 1/2, dndp = (pdp/n) cos 2 (3d (3, 



(37) 



where n is the principal quantum number. To the macroscopic accuracy that we are confining 
ourselves to, we can replace the total angular momentum j of the quasiparticle by its orbital 
angular momentum I, so that flj = I = psin (3. Taking also into account the fact that in our 
theory small values (3 and p ~ pf are important at low temperatures (pf — kfR, where kf is 
the wave number at the Fermi boundary) we can put 



~ df 

n 3 = Pf (3, g 3 = - 



9j, 



T = — 

dp 



T. 



(38) 
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In order to maintain the required macroscopic structure (see the Introduction) in this same 
limit, it is necessary, as was shown in [2], that 

9j = 9$. (39) 
The distribution function w s assumes the following form in the variables p and j3 (see Fig. 3): 



w s (p,(3) = { 



exp 



Pf-p 



exp I 7 I 



-l 



P~Pf<-ri 



-T7 



( 9Pj£ 

V 2 - 



< p - p f < 



ri , (40) 



exp 



(^) 



+ 1 



P-Pf>ri 



We also transform to the variables p and (3 in Eq. (25), and using the Jacobian dndl = 
(pf /ir)dpdp (see Eq. (37)) to replace the summation over the levels by an integration, we 
obtain 



S = 



oo poo 



oo JO 



2p f (3 [w s In w s + (1 - w s ) ln(l - w s )} ^dpd(5. 

71 



(41) 



Substituting (40) into (41) and integrating, we obtain the following expression for the entropy 
of one nucleon component of a spherical nucleus: 



S = 



7T 3 PfT 2 7T 3 



Pf 



4 g A g{de/dp\ f y 



t2 = 7T 3 R*m* 2 T2 
" 4 tfgpf ' 



(42) 



where m* = Pf/vf is the effective mass of a quasiparticle in the nucleus. Taking into account 
the fact that nuclear matter is a two-component system, we find 



b 
2 



S = ^T 2 , b 



7T 3 R A m* 2 
T H 4 



1 1 

+ 



Pf9N pjgz 



(43) 



Thus, in spherical nuclei (at low enough temperatures — see below) the entropy depends 
quadratically on the temperature, in contrast to the frequently assumed linear dependence. 
The latter is characteristic of the energy spectrum of a normal Fermi liquidl [9,10], but in nu- 
clear physics it can apparently be regarded as fairly reliably established only for nonspherical 
nuclei. 

The thermodynamic relations 



E = -T 3 , 



C =^ = bT 2 , S=—b 1 / 2 E 2 / 3 , 
dT ' 2 



(44) 



follow from Eqs. (43), so that the level density in a spherical nucleus depends as follows on 
the excitation energy: 



Pi.d.(E) oc exp[S(E)} = exp 



32/3 



2 b^E^j. 



(45) 
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Now let us calculate the moment of inertia of a spherical nucleus. To do this we shall 
examine the Racah-Mottelson states that have not only a definite seniority Sj, but also a fixed 
angular momentum projection Mj. Then the entropy of level j will be Sj = lnT(Mj, Sj, Qj), 
where the number of states is given by Eq. (24) (as before, we neglect the preexponential 
factor). The distribution functions for the quasiparticles and the angular momentum Mj are 
determined from the condition that the free energy F' = E — TS — EjN — ooM be minimum, 
where u is the angular velocity and M is the total angular momentum projection. Thus, the 
Mj dependent terms 

3M?T 



E 



n jSj (2n j - Sj 



uoMj 



(46) 



are added to the terms determined by formulas (25)-(27). Varying Mj yields the angular 
momentum projection 

Mj = J l ^ m > s ' ] (47) 
of the j level, so that the total angular momentum projection is given by 

M = ^E^( 1 -^)- ( 48 ) 

In this case the distribution function Wj depends on u, but for small uo we can neglect this 
dependence in Eq. (48) and use Eq. (40). The coefficient of uo in (48) is the moment of 
inertia: in the usual units (M in erg-sec and u in sec -1 ) it takes the form 

3 

Now we transform (49) to the variables p, (3 and substitute expression (40) for the distribution 
function; then after integrating over regions I — III we obtain 

/ = ^ ^ ^ T 2 (50) 
12 0* {de/dp\ f f ■ 1 ] 

In our calculation, however, we have not taken into account the effect of the second functional 
derivative of the energy with respect to the distribution function [9,10], which leads to an 
additional constant factor in Eq. (50). This factor can be found by calculating the moment of 
inertia of a normal Fermi liquid by a method similar to that used above. A simple calculation 
shows that the degree of degeneracy of a state with s fermions on level j having the angular 
momentum projection M differs from (24) only in the preexponential factor, which is not 
important for thermodynamic applications. This means that in the thermodynamic sense 
the seniorities s are fermion excitations with a dispersion law leading to distribution function 
(40). The moment of inertia of a Fermi liquid is therefore given by Eq. (49) with the ordinary 
Fermi distribution substituted for w s . 
Now we easily find 

/ - 4 h2p f 

FL 45n de/dp\/ 1 ; 
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Here, too, we have neglected the second functional derivative. Requiring that taking this 
into account should lead to the rigid-body value Ifl = (2/5)Nm n R 2 = (2 / 5) (2p 3 j \ / 97i~)m n R 2 , 
where m n is the nucleon mass, we find 



de 
dp 



h2p f (52) 



m 



R 2 ' 



i.e., de/dp\f should have the Fermi-gas value. Hence one of the factors de/dp\ f in Eq. (50) 
should also have the Fermi-gas value. Taking this into account, as well as the fact that 
nuclear matter is a two-component system, we obtain the expression 

_ vr 3 m* 2 m n R« ( 1 1 \ 

for the moment of inertia of a spherical nucleus. 

The function describing the distribution of the nuclear levels over the angular momentum 
J at fixed temperature has the form (see [11]) 

h 2 ~ ( h 2 J 2 \ ~ l 

wj = — Jexp — , J = J + -, (54) 

IT y 2IT J 2' K J 

so that the mean and mean square values of the angular momentum are 



7 HIT y 2 2IT 

Thus we have 

— 4 /-\ 2 

J 2 = -(j) , (56) 

and this relation can be used to check the correctness and completeness of the experimental 
data on level spins. We also note the equation 

« 9N9z{pf9N + Pf9z) 
which follows from Eqs. (55), (53), and (44). 



4 Discussion 

As we already mentioned, the formulas derived here are valid only at low temperatures. In 
order to find out how low the temperature must be, we take account of the fact that Eq. (39) 
was established only for angles satisfying the condition 

(58) 
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(see [2]). Hence the quantity 



! = \l— , (59) 
9Pf 

which is characteristic of our theory, marking, as it does, the point at which region II begins 
to have a finite extent along the e or p axis (see Figs. 1 and 3 and Eq. (40)), should satisfy 
condition (58). This gives the criterion 

r < § or T < fc, (60) 
8 

where fe = de/dp\j g/8 is the characteristic width of the diffuse region of the Fermi distri- 
bution caused by the residual reaction (see [2]). The values of 5e found in [2] range from 
1.9 MeV (heavy nuclei) to 5.2 MeV (light nuclei) for neutron magic nuclei and from 1.6 MeV 
(heavy nuclei) to 3.8 MeV (light nuclei) for proton magic nuclei. However, it is difficult at 
present to judge the part played by another criterion, 

T « T c , (61) 

implied by the theory, where Tq is the Curie temperature at which the spherical nucleus 
passes from the magic phase to the normal phase, i.e., becomes nonspherical (see [3]). We 
can only assert that condition (61) becomes more important, and may become dominant, 
for spherical nuclei that lie close as regards number of nucleons to the phase transformation 
point. The 7-ray spectra and radiative widths of nuclei that are spherical in the ground 
state confirm this (see [5,12]). 

We also note that all the equations we have derived relate the average values of the 
thermodynamic quantities to one another, i.e., they are valid only when the fluctuations are 
small. However, the temperature fluctuations [10] 

AT ^7I (62) 

amount to ~ 0.6 MeV for 208 Pb, so that the thermodynamic formulas can only serve as 
rough estimates at the excitation energies now being investigated for spherical nuclei. For 

example, it follows from Eq. (57) that J 2 ~ 4 for 208 Pb at an excitation energy of 5 MeV, 
whereas Eqs. (55) with the rigid-body moment of inertia and temperature T ~ 0.4 MeV 

(which is apparently characteristic of nonspherical nuclides; see [12]) give J 2 ~ 160. As well 
as can be judged from the available data [13] on levels of 208 Pb, however, the experimental 

value is J 2 ~ 10, which is in qualitative agreement with the theory under discussion. 

Although there is no thermal phase transformation in the present theory, the macroscopic 
structure of the residual interaction as a function of the orbital quantum number I together 
with the fact (a macroscopic quantum effect) that the spherical nuclear configuration is 
unstable in the absence of interaction between the quasiparticles (see [6]) enables us to 
suggest a possible reason for such a transformation. In fact, as the temperature rises region 
II shifts to the right on the / axis and more and more quasiparticles fall into regions I and III, 
where they conform to the ordinary Fermi distribution. Hence a phase transformation from 
the magic phase to the normal phase should take place when there are no longer enough 
quasiparticles left in region II to stabilize the spherical shape of the nucleus. 

I express my deep gratitude to V.G. Nosov for numerous discussions and valuable remarks. 
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